Sequential Tunneling vs. Electron Correlation in Multiple Photo-ionization 
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We take advantage of the information provided by use of elliptical polarization in a recent two- 
electron release time experiment [A.N. Pfeiffer et ai, Nature Physics 7, 428 (2011)]. This allows 
a comparative test of the currently dominant conjectures regarding independent-electron tunneling 
theory vs. fully electron-correlated classical release theory to describe electron release in strong-field 
double ionization. 

PACS numbers: 32.80.Rm, 32.60,+i 



Atoms exposed to very strong laser fields are sub- 
ject to ionization. With some experimental support [l|, 
quantum tunneling is generally conjectured as the phys- 
ical process that is responsible for such ionization: an 
electron can be freed by "tunneling through" the poten- 
tial barrier generated by the combination of the atomic 
Coulomb potential and the laser electric potential. 

Theories have been proposed to predict the tunneling 
rate of an atom exposed to a laser field 0^3 • These the- 
ories are based explicitly on an assumption that the laser 
electric field is weaker than the over-barrier field of the 
ionizing electron. The errors of these theories are diffi- 
cult to estimate, but certainly grow substantially as the 
laser field becomes stronger, while the theories continue 
to be used for guidance. However, their applicability is 
now beginning to be open to quantitative test by compar- 
ison with recent experiments that record electron release 
times using some of the advantages provided by elliptical 
polarization of the laser field (for an overview, see [5|). 

To extend the applicability of the tunneling theories 
and at the same time keep the simplicity and convenience 
of the analytical form of the well-known Ammosoy- 
Delone-Krainov tunneling formula Tong and Lin Q 
proposed an augmented quantum tunneling (AQT) for- 
mula in which an additional exponential correction factor 
is added to account for the strong-field saturation effect. 
The correction factor in the AQT formula contains an 
empirical parameter a to be determined for each atomic 
and ionic species. 

Under the usual tunneling conjecture, an electron's re- 
lease time is not really open to question. Because of 
the exponential dependence of the emission rate on laser 
field strength predicted by tunneling theories, the elec- 
tron is almost certainly released at the peak of the laser 
pulse. Modern laser technologies can now routinely gen- 
erate laser fields that are strong enough to ionize multiple 
electrons from an atom. These laser fields can be near to 
or higher than the over-barrier field of the first electron or 
of the second electron, so an electron can be emitted ear- 
lier than the peak of the pulse, and determining multiple 



FIG. 1: Comparison of the AQT predictions (solid lines) of 
emission times with experiment (dots with error bars). The 
left panel is for the 7fs pulse and the right panel is for the 
33fs pulse. Blue color is for first ionization and green color is 
for second ionization. Adapted from [7|. 



emission times becomes a viable route for testing. 

The AQT formula has been used by Pfeiffer et al. Q 
in connection with their measurements of the sequential 
release times in double ionization of Ar, using ellipti- 
cal polarization. To deal with two-electron ionization, 
the AQT theory uses a single-active-electron (SAE) ap- 
proximation, which assumes that only one electron is ac- 
tive at each time, with other electrons playing no role 
but screening the nucleus. Only after the first electron 
is ionized does a second electron become active. Com- 
parison of the AQT predictions with experiments show 
good agreement with the observed release time of the 
first electron, while they predict substantially later times 
than observed for the second electron. This is clearly 
seen in Fig. [1] It has therefore been conjectured that 
the SAE (or independent-electron) approximation breaks 
down, meaning that electron-electron correlations play a 
role. However the specific role that they play has not 
been explained. 

The same conjecture of required e-e correlation has 
also been advanced for another experiment, with a simi- 
lar reasoning Q. In that case, the tunneling theory fails 
to explain an oscillation feature of the measured ratio 
between parallel and anti-parallel sequential double ion- 
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ization (SDI) emissions (i.e., between SDI with the two 
electrons released in the same or opposite directions). 
Interestingly, a fully correlated classical ensemble simu- 
lation based on solution of the relevant time-dependent 
Newton equations (TDNE) qualitatively reproduces the 
oscillation feature [9j . The physical role of the e-e corre- 
lations was again not clearly identified. 

Moreover, Zhou et al. [lOj recently used a classical 
ensemble simulation, again employing TDNE solutions 
of the relevant equations, and reported the first success- 
ful reproduction of the measured electron release times. 
They modelled each atom in their ensemble as contain- 
ing two fully correlated electrons able to interact steadily 
with each other, as well as with the laser field and the ion. 
But the specific role of the e-e correlations leading to the 
experimental agreement were not identified or analyzed. 

From this short review, one sees that conjectures have 
been advanced that e-e correlations must necessarily be 
present on either of two different grounds: because (1) 
the independent-electron SAE-tunneling AQT approach 
disagrees with the experiments, and/or (2) the fully- 
correlated classical simulation agrees with the experi- 
ments. 

Of course, making explanatory conjectures such as the 
need for e-e correlation should be encouraged. But the 
lack of an understanding of the role for such correlations 
substantially weakens the conjectures. There is a reason 
for this. As we report here, systematic calculations show, 
surprisingly, that the need for e-e correlation on the one 
hand, or quantum tunneling on the other hand, can both 
be discarded in treating electron release timing. With- 
out either one, quantitative agreement with experimental 
release times can still be achieved. 

Classical ensemble calculations (traditionally fully- 
correlated) have previously been described in detail 
[ill [l2j]. Recently the classical ensemble method has 
been extended to include elliptical polarization 13- ljj 
and good agreement with experiment has been achieved 
0, Ea]. In the calculations reported here, the same 
method is implemented, however while enforcing SAE 
evolution. This is done by numerically permitting only 
one electron to be actively involved in the ionization pro- 
cess, and then after the first electron is removed by the 
laser field a second electron is allowed to become active 
and subject to removal from the residual ion. This pro- 
cess can be extended to additional electrons as desired. 

Without difficulty, the Newtonian equations of motion 
can be integrated numerically under the SAE constraint 
and the positions and the momenta of the two electrons 
are recorded at each time step until the end of the pulse. 
At the end of the pulse, ionization trajectories with both 
electrons at least 10 a.u. away from the ion core can be 
collected and further analyzed. Typical two-electron re- 
lease trajectories are shown in Fig. [2j corresponding to 
a time delay of an integer number of cycles between the 
two emissions (left) and of an odd multiple of half cycles 
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FIG. 2: Two typical SDI trajectories with elliptical polariza- 
tion. The left panel shows a time delay of an integer number 
of cycles between the two emissions, and the right panel shows 
a time delay of an odd multiple of half cycles between the two 
emissions. 



between the two emissions (right). Trajectories as shown 
in the left and right panels usually lead, respectively, to 
larger and smaller net ion momenta along the minor di- 
rection, analogous to the non-Z and Z trajectories found 
with linear polarization [17]] . Detailed explanations can 
be found in [13j . 

Some specific details of the classical SAE simulation 
procedure are described as follows. Initially an ensem- 
ble of classically modeled atoms with only one electron 
is generated. The energy of this electron is set to the 
negative of the first ionization potential (noted as I Pi ) 
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(1) 



where pi and i"i are the momentum and the position of 
the electron. Note that the ion core Coulomb potential 
has been softened to avoid numerical singularities [l2j ]. 
Argon, which is the target atom used in the Pfeiffer ex- 
periment, has I Pi = 0.58 a.u. Given the energy, then the 
position and the momentum of the electron are randomly 
assigned. Therefore the initial ensemble is a collection of 
model atoms with the same energy but different position 
and momentum configurations. 

Then a laser pulse is turned on. We use the same laser 
parameters as in [tJ- The laser field is elliptically polar- 
ized with ellipticity value e = 0.77. The pulse envelope 
has a Gaussian shape of duration (FWHM) 7fs or 33fs. 
The laser field has the form 



E(t) = 



/(*) 



[e x sin(wi + ipo) + e y e cos(wi + <p )] 



(2) 

where f(t) = E exp (— < 2 /(2er 2 )) is the pulse envelope. 
Eq, cj, and e are the amplitude, angular frequency, and 
ellipticity of the field. The carrier envelope phase is noted 
as ipo, which is not stabilized and varies from pulse to 
pulse. Consistent with Refs. [HO, [HI], the x direction 
has been selected as the major polarization direction and 
the y direction as the minor polarization direction. 
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FIG. 3: Left: a correlation spectrum of pi r and p2 r for a 33fs 
pulse with peak intensity 4.5 PW/cm 2 . Right: reduced dis- 
tributions of pir and P2r, under the same condition. Heights 
of peaks have been normalized to unity. 



The Hamiltonian of this electron is 



H = H(t) 



Pi 



2 y/vT+T 



•ri-E(t). 



(3) 



If the laser field is strong enough, the electron can be 
ionized and driven away from the vicinity of the ion core. 
Only then, in the SAE a ppr oximation, is a second elec- 
tron allowed to be active |18| . In our approach, a second 
electron is "created" near the ion core when the first elec- 
tron reaches a distance of 10 a.u. A slight different choice 
of this distance will not affect our discussions at all. 

The initialization of the second electron is similar to 
that of the first electron, and the energy of the second 
electron is 



2 >/rf+l 



= -IPi 



(4) 



where IP% — 1.02 a.u. for the second electron. Given 
the energy value, the position and the momentum of this 
electron are randomly assigned. 

The second electron reacts immediately to all the forces 
on it. These forces include the Coulomb attraction from 
the ion core, the laser electric force, and the repulsive 
force from the first electron, which is already negligibly 
weak. The motion of the second electron is also governed 
by the classical TDNE. If the laser field is strong enough, 
the second electron can be ionized. 

The momenta of the two ionized electrons at the end of 
the pulse are recorded as {pi x ,Pi y ,Piz} and {p2x,P2 V ,P2z}- 
Two new scaled momenta are defined Q: 



Pi 



J(l + e*)(pl x /e*+pl y ); (5) 

P2r = yj(l + e*)(? 2x /e*+p> y ). (6) 

A correlated spectrum of p\ r and p2r is shown in the left 
panel of Fig. [3] for the 33fs pulse with peak intensity 4.5 
PW/cm 2 . It is diagonally symmetric due to the fact that 
no information about the ionization order is pre-known 
experimentally. Nevertheless, because the first ionized 



electron has a smaller scaled momentum p r than the sec- 
ond ionized electron, the diagonally symmetric correla- 
tion spectrum can be separated into individual distribu- 
tions of the first and second ionized electrons, as shown 
in the right panel of Fig. [3] The peak positions of the two 
distributions will be used to calculate the two emission 
times. 

Given that both ionizations happen before the peak of 
the pulse, the two emission times can be uniquely deter- 
mined from the peak positions of the distributions of p\ r 
and p2r (treated as Gaussian) [lj| : 



t; = —o\ / 2 In 



Eq 
UJPr 



for i = 1,2. 



(7) 



The results of the calculated ionization times as a func- 
tion of intensity are shown in Fig. [4j for both the 7fs and 
the 33fs pulses. Focal volume effects have been taken 
into account [2(|, as done in Q and (loj . Quantitative 
agreement is highly satisfactory over the full experimen- 
tal intensity range. 

The physical reason for the ability of the classical sim- 
ulation to match experimental results, as shown in Fig. 
21 while the AQT theory cannot, as shown in Fig. [TJ re- 
mains an important open question. The fact is that both 
are based on the SAE approximation, so their key dif- 
ference is simply rejection or adoption of initiation by 
tunneling. The other interesting open issue is in the 
agreement of our results not only with the experiments, 
but also with the calculations of Zhou et al. Those fully 
include e-e correlation, which our results suggest are sur- 
prisingly unimportant. 

We need to emphasize that the above argument by no 
means proves that the classical model is correct while 
the AQT theory is wrong. One may equally argue that 
the AQT theory may be correct, while the disagreement 
between the Pfeiffer experiment and the AQT predic- 
tion may originate from e-e correlations or multiclec- 
tron effects excluded from the AQT theory. We cannot 
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FIG. 4: Comparison of TDNE predictions of emission times 
with experiment, obtained with new calculations, as explained 
in the text. The left panel is for the 7fs pulse and the right 
panel is for the 33fs pulse. Dots with error bars are the exper- 
imental data adapted from [7J. Lower and upper solid lines 
are the classical SAE predictions of the 1st and the 2nd ion- 
ization, respectively. 
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deny this possibility but we do not like it. First, one 
always prefers simpler theories without many complica- 
tions. And second, the conjectured e-e correlations or 
multiclcctron effects arc quite difficult to identify, partic- 
ularly as the results of the fully correlated theory of Zhou 
et al. is entirely in accord with our completely uncorre- 
cted results. 

To summarize, in our calculations we have discarded 
any remnant of either of the conjectures (1) or (2) iden- 
tified at the outset. This was done quite simply. We 
adopted the usual classical ensemble theory, in order to 
strictly preclude tunneling. And then we numerically 
prevented e-e correlation. In this way we simultaneously 
tested both the importance of tunneling and the neces- 
sity of electron correlation. Classical calculations have 
regularly provided useful trajectory-based insights into 
strong-field ionization dynamics 17J, l2l| - |25j , and we have 
no reason to suspect anything different here. 

In conclusion, we have shown that the Pfeiffer et al. 
experiment on electron release times of SDI can be quan- 
titatively reproduced by a tunneling-free, correlation-free 
classical analysis. The physical reason for the difference 
between the two SAE theories, our classical simulation 
and the AQT result, is yet an open question to answer. 
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